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Abstract: We develop a consistent semiclassical description of field-theoretic collision-indnced 
tnnneling at arbitrary high collision energies. As a playgronnd we consider a (l + l)-dimensional 
false vacnnm decay initiated by a collision of N particles at energy E, paying special attention 
to the realistic case of = 2 particles. We demonstrate that the cross section of this process 
is exponentially snppressed at all energies. Moreover, the respective snppressesion exponent 
F]sf{E) exhibits a specific behavior which is significant for onr semiclassical method and as- 
snmed to be general: it decreases with energy, reaches absolnte minimnm F = Fmin{N) at a 
certain threshold energy E = Ert{N), and stays constant at higher energies. We show that 
the minimal snppression Fmin{N) and threshold energy can be evalnated nsing a special class 
of semiclassical solntions which describe exponentially snppressed transitions bnt nevertheless 
evolve in real time. Importantly, we argne that the cross section at energies above Ert{N) is 
compnted pertnrbatively in the backgronnd of the latter solntions, and the terms of this per¬ 
turbative expansion stay bounded in the infinite-energy limit. Transitions in the high-energy 
regime proceed via emission of many soft quanta with total energy Ert] the energy excess E — Ert 
remains in the colliding particles till the end of the process. 
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1 Introduction 

Exponential suppression of probabilities precludes direct observation of extraordinary tunneling 
phenomena such as baryon number violation in instanton-like electroweak transitions [1-3] 
or spontaneous decay of allegedly false Higgs vacuum [4-8]. Quantum mechanical intuition 
suggests, however, that tunneling probabilities grow with energy. Indeed, tunneling phenomena 
of the above sort occur at higher rates [9, 10] in two-particle collisions: 

a{E) oc , (1.1) 

where the suppression exponent^ F{E) decreases with collision energy E, while g is a. small 
coupling constant. The central question is whether the exponential suppression in Eq. (1.1) 
disappears at sufficiently high energies and, if it does not, what is the value of the suppression 
exponent at E —)■ +cx). 

This question is surprisingly nontrivial. Field-theoretic tunneling involves barriers of hnite 
heights Ecb given by the energies of the critical bubble [12-15] and sphaleron [16, 17] in scalar 
and gauge theories, see Fig. la. Nevertheless, the respective collision-induced transitions cannot 

^Dubbed “holy grail function” [11]. 
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Figure 1. (a) A sketch of the potential barrier Epot[4)] between the vacua (j) = 0 and cj) = 

(b) Scalar potential V{p). (c) Suppression exponent Fjv(-E) at high energies. 


become unsuppressed^ at > Ecb [20, 21], Consider e.g. massless fermions -0 and -0 conpled 
with small Yukawa constant Y to the scalar sector of the theory.^ Their contribution to the 
scalar self-energy n(Q^) obeys dispersion relation [21], 



where atot{E) is the total annihilation cross section '0-0 —)■ anything at the center-of-mass en¬ 
ergy E and we ignore irrelevant scalar masses. If the fermion annihilation leads to unsuppressed 
over-barrier transitions ai E > Ecb, the self-energy n(Q^) in Eq. (1.2) receives large nonper- 
turbative contributions at small Q^. This would contradict the standard perturbation theory 
which is valid at low energies. Thus, the cross section (1.1) of collision-induced tunneling is 
exponentially suppressed at E > Ecb- 

In this paper we describe collision-induced tunneling at high energies. To be specihc, we 
consider false vacuum decay in the (1 -|- l)-dimensional scalar held model 



VW 


(1.3) 


with false and true vacua at 0 = 0 and (p = <p+, respectively; the scalar potential V{(p) is 
shown in Fig. lb. We work at weak coupling, -C 1. At zero energy the vacuum 0 = 0 
decays spontaneously via formation of an expanding bubble with (p ^ (p^ inside. Below we 
study the same decay accompanied by a collision of N 0-quanta at energy E. We compute the 
suppression exponent En{E) of the corresponding inclusive cross section at high energies. 

Our numerical result‘d for E]sf{E) is presented in Fig. Ic. This function decreases with 
energy, reaches minimum Emin{N) at E = Ert{N) and stays constant at higher energies. With 
the aid of the Rubakov-Son-Tinyakov conjecture [22] we extrapolate results to the two-particle 

^Gravitational interactions are completely different in this respect [18, 19]. 

^It is straightforward to generalize this argument to other setups. 

'^In numerical calculations we use a specific form of V{p) which is not significant at the moment. 
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initial state and find similar behavior of the exponent F{E) = F 2 {E) in Eq. (1.1) (upper panel 
in Fig. Ic). We conclude that collision-induced false vacuum decay is exponentially suppressed 
at arbitrary high energies. 

Energy-independent suppression exponent F = Emin of collision-induced tunneling at high 
energies was proposed in Ref. [23] and observed in toy models of Refs. [24, 25]. We hnd the 
same behavior in the full-fledged held-theoretic model. 

Besides, we demonstrate that induced false vacuum decay at the threshold E = Ert{N), 
despite being exponentially suppressed, is described by one-parametric family of complex semi- 
classical solutions 0 = 4>rtix) evolving in real time t = . These solutions were introduced in 

Ref. [26] under the name real-time instantons. They satisfy complexihed classical held equa¬ 
tions with certain boundary conditions in the asymptotic past and future. The initial particle 
number N parametrizing the solutions enters the conditions at t —)■ — oo. One can show [26] that 
the minimal suppression Fmin{N) and respective energy Ert{N) are computed as functionals 
on the real-time instantons. 

We argue on general grounds that if the real-time instantons exist for some collision- 
induced process, the respective suppression exponent is energy-independent and equal to 
Fmin{N) aX E > Ert{N). Thus, Ending the family of these solutions in a given model, one 
obtains the exponent F]\f{E) in the entire high-energy region E > Ert{N). 

Most importantly, we demonstrate that the real-time instantons serve as backgrounds for 
the long-awaited [23, 27] perturbative description of collision-induced tunneling at high en¬ 
ergies, and the respective corrections are bounded in the high-energy limit. This remarkable 
feature is in sharp contrast to the properties of perturbative expansions in Euclidean back¬ 
grounds [9, 10] which blow up [11, 28] [3] at E ~ E^b- Let us explain the difference by considering 
scatterings of particles at energy AE in different backgrounds. At a crude level the quantum 
particles can be regarded as small-amplitude high-frequency waves 6(j) oc added to the 

background. In the Euclidean case 6(() grows as with r = —it, and nonlinear backreaction 
effects become essential at high AE. In other words, the perturbative expansion in 50 breaks 
down. In the opposite case of the real-time instanton the waves 50 evolve adiabatically and do 
not change the soft background. Scattering of these waves can be described perturbatively. 

Using the above observation, we propose a working perturbative scheme [29] for evaluating 
the inclusive cross section of collision-induced tunneling as series in in the most interesting 
case of high energies E > 77^(2) and two initial particles. The receipt is as follows. One starts 
by developing a formal perturbation theory in the background of a real-time instanton (j)rt{x) 
with parameter N = No. Namely, one considers the {n -|- 2)-point Green’s function 

Grt = {'^rt\4>{xi) . . .(j){Xn+2)\'^o) = j "^0 4^^ [0] 0(Xi) . . . 0(a;„+2) Tq [0] (1.4) 

between the false vacuum Tq and the most probable hnal state^ of this real-time instanton. 
One substitutes 

(p{x) = 0rt(a;) +g5(j){x) 

coherent state to be specified in the main body of the paper. 


3 



into Eq. (1.4) and evaluates the path integral over fluctuations Sep as series in g. The terms in 
these series depend on the auxiliary parameter Nq of the real-time instanton which characterizes 
the background configuration. Using the LSZ reduction formula, one extracts from the Green’s 
function the amplitude of the process 2 —>■ n + where the final state contains n particles 
on top of Trt. We argue that the final states of this kind form a complete set. Summing over 
them, one finds the inclusive cross section of the collision-induced process. The final result is 
obtained in the limit Nq —)■ 0. 

Our method is different from the standard instanton perturbation theory [9, 10] in several 
respects. First, we use the real-time instanton as a background. This guarantees stability of 
the perturbative expansion at high energies. Second, we introduce an additional parameter 
of (prt and send Aq —0 in the end of calculations. Indeed, the real-time instanton with 
Nq = 0, as we argue in the main body of the paper, is the saddle-point configuration for the 
path integral (1.4). Formally it describes transition from vacuum: two initial particles of the 
process are represented by two 0-factors in the integrand of Eq. (1.4) which do not affect the 
saddle-point solution. On the other hand, (prt{t, x) turns out to be singular at Nq = 0. We 
therefore work with smooth saddle-point configurations at Aq > 0 and recover® the correct 
results in the limit Aq —?• 0. Third and Anally, although below we evaluate only the leading 
inclusive suppression exponent F = Fmin at A > our perturbative approach can be used 
for the prefactor and exclusive cross sections at high energies. 

With the help of the perturbative method we identify^ the 
dominant mechanism of the collision-induced transitions at A > 

Art(2) and A = 2, see Fig. 2. We observe that the colliding par¬ 
ticles emit many soft quanta which form a bubble of true vacuum 
with energy Art(2); the energy excess A — Ert{2) remains in the 
initial particles till the end of the process. 

From a general prospective our results support exponential 2 Transitions at 

suppression of collision-induced tunneling at arbitrary high ener- ^ ^ Ert{‘2) 
gies and per se put on shaky ground proposed searches for non- 

perturbative phenomena at future colliders [32] or in cosmic ray events [33, 34], cf. Ref. [35]. 
For example, it was found in Refs. [36, 37] that the suppression exponent of electroweak baryon 
number violation in two-particle collisions is almost energy-independent at A ~ 15 TeV. If the 
minimum A = Ert{‘2) is somewhere near this point, the respective cross section is suppressed 
at all energies by a deadly factor TeV)/aw ^ where a-w is the electroweak cou¬ 

pling and we took numerics from Refs. [36, 37]. To get reasonable probabilities, one should 
consider models with tunneling rates raised by dynamical mechanisms [38-41], resort to strong 
coupling [42, 43] or exotica [44]. 

We argued in Ref. [45] that the collision-induced false vacuum decay in (1 + 1) dimensions 
turns into production of kink-like soliton pairs from particles once the energy densities of the 

®This makes our procedure a regularized version of Landau method [30] [24, 31] of singular semiclassical 
solutions. 

^For simplicity we consider the process in a large finite volume. 
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two vacua are leveled, V{0) — V (0+) —)■ 0. One expects that the properties indicated above hold 
in this limit. In particular, the cross section of creating a pair of solitons from two particles 
is exponentially small at all energies, and the suppression exponent of this process does not 
depend on energy above a certain threshold. 

This paper is organized as follows. In Sec. 2 we recall perturbative expansion in the 
background of a Euclidean bounce. This technique reproduces exponentially growing collision- 
induced cross section at low energies but breaks down at E > Ecb- We proceed with moderate 
energies in Sec. 3 and demonstrate how the bounce at E ^ Ecb is connected with the real¬ 
time instantons at E = Ert{N) > Ecb- The latter semiclassical solutions and transitions at 
E > Ert{N) are considered in Secs. 4 and 5, respectively. In Sec. 6 we show that the same term 
of perturbative expansion which led to exponential growth of the cross section with energy 
in perturbation theory around the bounce, gives subdominant and exponentially decreasing 
contribution in the background of the real-time instanton. We propose a working perturbative 
description of collision-induced tunneling at N = 2 and high energies in Sec. 7. Our results are 
summarized in Sec. 8. 


2 Perturbative expansion in the background of a bounce 


Let us explain the difficulties with the collision-induced tunneling by reviewing its low-energy 
description [9, 10] in the model (1.3). We will also introduce terminology and sharpen contrast 
with the high-energy transitions. 


Spontaneous decay of false vacuum at E = 0 is described 
by the celebrated bounce solution [14, 15] (j)b{x), see Fig. 3. The 
latter has Euclidean and Minkowski parts representing nucleation 
of a true vacuum bubble and its expansion to inhnite size. Note 
that the bounce is Lorentz-invariant i.e. depends on = x^x^. 
At large negative x^ it satishes the Klein-Gordon equation in the 
false vacuum and therefore behaves as 

0fe(x) —)■ ^ Ko{m\/—x‘^) as x^ —)• —oo . (2.1) 

27r 

Here m is the mass in the false vacuum 0 = 0. In Eq. (2.1) and 
below we exploit universal complex time f = x° which is real and 
pure imaginary in the respective parts of the contour in Fig. 3. 
Parameter Cb is related to the bubble size Rb- cl)b is of order 1 at 
X ~ Rb- In the thin-wall limit mRb ^ 1 one obtains Cb oc 
where the asymptotics of the Bessel function in Eq. (2.1) was 
used. 

To compute the amplitude of collision-induced tunneling, we 


l\ 




Figure 3. Bounce 0b(x). 
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consider the Green’s fnnction 


Qb = (d'b|0(xi)... 0(a;n+2)|d'o) 

= 1^0 n[^] 0(xi)... 0(x„+2) ^o[0] , (2.2) 

where |\l/o) and |\l/b) are the false vacnnm and the final state of its decay at E = 0, respectively. 
The latter describes expanding bnbbles. 

The idea [9, 10] is to evalnate the path integral in Eq. (2.2) in the saddle-point approxi¬ 
mation treating ((){xi)... 0(a;„+2) as a prefactor. The relevant saddle-point confignration is the 
bonnce (pb- In the leading order one obtains, 

Gb = Ab j (fxo (j)b{xi - xo)... 0fe(x„+2 - xo ), (2.3) 

where we introdnced the bonnce position xo, ignored irrelevant saddle-point determinant and 
denoted the bonnce amplitnde by Ab = tho[06]. Note that the action S and zero- 

energy states To, are translation-invariant® and therefore independent of Xq. The bonnce 
amplitnde is compnted in Refs. [14, 15]: with exponential precision 

|^^|2 = , (2.4) 

in particnlar, |Tft[0b]|, |To[06]| ~ 1. We remark that corrections to Eq. (2.3) can be evalnated 
pertnrbatively nsing the bonnce as a backgronnd: one snbstitntes (f) = (j)b + gS(j) into the path 
integral and calcnlates it as series in 

Importantly, the large-x asymptotics of the bonnce in Eq. (2.1) reprodnces correct residnal 
at k‘^ = m?‘ of its Fonrier transform, 

/ XCu 

d^xe*^'’’’06(x - Xo) = Yo - 9 -^ + regnlar part , (2.5) 

+ le 

where we recalled that Kq in Eq. (2.1) is the Feynman propagator in two dimensions; we denote 
k-x = k^x^ and assnme e —)• -|-0. The residue Cb in Eq. (2.5) is fc-independent or “point-like” [3]. 
This property is specific to Euclidean solutions, as we argue below. One immediately obtains 
the 2 —)■ n transition amplitude from the LSZ formula, 

/ \ n+2 

A2^n = Abl^] . (2.6) 

To derive this expression, we Fourier-transformed Eq. (2.3), extracted the on-shell residues 
(2.5) and collapsed the integral over xo into the ^-function® representing the energy-momentum 
conservation. Factors g in Eq. (2.6) compensate for non-canonical normalization of kinetic term 
in Eq. (1.3). 

®Here and below we assume asymptotic limits tij —>■ =Foo when all quantities oscillating with the initial ti 
or final tf times can be dropped. 

®It is absorbed in the phase space volume, as usual. 


6 




The amplitude (2.6) exponentially grows with energy E if the most probable hnal state 
contains n ^ E/m nonrelativistic particles. To confirm this guess about the most probable 
state, we derive the inclusive cross section in Appendix A, 


a{E) = ^ / |-4,2^n| dUn = \Ab 


2 \ 


d^Xe 


27vg 


(2.7) 


where the prefactors are ignored, n„ is the n-particle phase space volume and P'^ = {E, 0) is 
the total momentum in the center-of-mass frame. The variable in Eq. (2.7) can be regarded 
as a typical Compton wavelength of the hnal particles: the latter become nonrelativistic at 
A 3> m~^. At g ^ 1 the integral in Eq. (2.7) is evaluated in the saddle-point approximation. 
The extremum of the exponent is achieved at 


A^ = {-2zT, 0) , 


T 


J_j f A 

2m \2g‘^EVw) ’ 


( 2 . 8 ) 


where the asymptotics of the Bessel function was used. This gives 


criE) = b®-2imS[0i,] 


(2.9) 


see Eq. (2.4). In the thin-wall limit V{(j)+) —)■ 0 we substitute Cb oc and obtain T ^ Rb 
with corrections proportional to log(mi?b). The cross section (2.9) in this limit coincides with 
that in Refs. [46, 47], see also Refs. [48-50]. 

The result (2.9) demonstrates exponential growth of the cross section with energy. It 
involves nonrelativistic hnal particles when A^ ~ T is large i.e. at E where 

the exponent comes from Cb- Moreover, the relativistic regime is never reached because at 
E > m/g"^ weak coupling expansion becomes unreliable. Indeed, relative corrections to the 
leading-order result (2.3) are estimated^*^ [3] as g'^n^, where v? comes from combinatorics; they 
are already large at n ~ E/m ~ One concludes that correct description of collision- 

induced tunneling at E > Ecb ~ m/g"^ should incorporate backreaction of the hnal-state 
particles on the semiclassical solution; we will pursue this approach in the next Section. 

Let us point at two specihc features of the low-energy calculation. First, the bounce residue 
Cb in Eq. (2.5) does not depend on k as k ^ +(X). Second, the hnal state '^b has zero energy 
and factorizes in Eq. (2.3). We will see that these two properties do not hold for perturbative 
expansion about the relevant high-energy solutions. 


3 Prom Euclidean to real-time solutions 

We demonstrated that collision-induced transitions are no longer described by the bounce at 
E > Ecb- Since our interest lies in high energies, we set this background aside and search 

Sophisticated resummation [51, 52] shows that the true expansion parameter is g^n 1. 
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for true semiclassical solutions describing the false vacuum decay in the iV-particle collisions. 
Consider the inclusive cross section, 


ME)= Y1 E, N) 


^-FN(E)/g^ ^ 


(3.1) 


where U is the evolution operator, tij —)■ =Foo, we ignored the initial flux in the prefactor and 
introduced the suppression exponent F]y{E) in the approximate equality. The sum in Eq. (3.1) 
runs over all initial states Tj with energy E and multiplicity N in the false vacuum and hnal 
states Tj containing a bubble of true vacuum. Importantly, a^iE) coincides at N = 2 with 
the two-particle collision-induced cross section o'(E) and can be computed semiclassically at 
3> 1. Moreover, Rubakov-Son-Tinyakov conjecture [22] states that the suppression exponent 
En {E) does not depend on at A^ -C see Refs. [53-56] for conhrmations. This means 
that the two-particle exponent E{E) in Eq. (1.1) is obtained by extrapolating the semiclassical 
result for En{E) to g'^N —)■ 0. 

In direct semiclassical approach one writes a path integral for (Jn{E) and evaluates it at 
g ^ 1 using the saddle-point conhguration 0s(^, In general, this configuration is complex. 
The saddle-point equations for 0^ are derived in Ref. [22], see the summary in Fig. 4a. Like 
the bounce, this conhguration satishes the classical held equations^^ 55/(50 = 0 or 

(5/ - 9^)0, = -R'(0,) (3.2a) 


along the complex time contour ABCD in Fig. 4a. It contains an expanding bubble in the 
asymptotic future where the solution is real. 


Im0s, lmdt4>s —t 0 as t —)■ +cx) . 


(3.2b) 


Peculiarities of 4>s{x) are related to the nontrivial initial state in Eq. (3.1). Along the part AB 
of the contour the saddle-point solution describes motion of the particles prior to the collision: 
in the asymptotic past it reduces to free waves 


0s 



as t —)■ —oo , 


(3.2c) 


where = {uk, k), ul = and the limit is taken along the time contour of Fig. 4a. 

In Eq. (3.2c) we introduced the classical counterparts a^, bl of the annihilation and creation 
operators with relativistic normalization. They are related by the initial condition 




^-2u>kT-e^^ 


(3.2d) 


involving two Lagrange multipliers T and 9 due to hxation of energy E and the number N of 
colliding particles. The latter quantities are given by the standard expressions 

= = (3.3) 

the spatial boundary we impose the standard energy-conserving condition dx4>s —>■ 0 at a; —> ±cx). 
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Figure 4. (a) Contour in complex time for x). Saddle-point equations and boundary conditions 
at t —)• ±oo are written near the respective parts of the contour. The nearest singularity t* of the 
solution is marked by the crossed circle attached to the branch cut (thin double line), (b) Physical 
quantities in the model (3.6). (c) Suppression exponent at low energies. 


In the limit T —)■ +oo Eq. (3.2d) reduces to vacuum condition ak = 0 which corresponds to 
E = N = 0. In this case 0s(x) coincides with the bounce solution (j)b{x)- At hnite T and 6 
the saddle-point solution describes transition at nonzero E and N. In what follows we solve 
equations (3.2) and relate (T, 9) to (E, N) by Eq. (3.3). 

Given the saddle-point conhguration (j)s{x), one evaluates the suppression exponent [22] 


E^iE) = /(2Im^[0s] - 2ET - NO) + Im 


dx (j)sdt(j)s 


t=ti 


(3.4) 


where the last three terms are the initial-state contributions. Importantly, the method of 
Lagrange multipliers implies Legendre transform 


dEFn(E) =-2gX , dMFn(E) =-g'^B , 


(3.5) 


which demonstrates that T and 9 are proportional to the derivatives of Ej^{E). 

We solve the boundary value problem (3.2) for different T and 9 numerically. To this end 
we specify the scalar potential in dimensionless units, 


1 /( 0 ) 


£ 

2 


l-vW 



(3.6) 


where W{x) = {x + + x^), u = 0.4, and the energy density 1/(0+) = —0.4 of the true 
vacuum is set by tuning v ~ 0.84. Function (3.6) is plotted in Fig. lb. It is almost quadratic 
at 0 < 0+/2 and nontrivial at larger 0, so that waves in Eq. (3.2c) remain linear^^ almost up 
to their collision point. 

We numerically computed the physical quantities for the potential (3.6), see Fig. 4b. To 
this end we have found the bounce 0b(t, x) and the critical bubble 0cb(x), see Ref. [15] for 

^^Long nonlinear evolution in other models is costly for numerical computations. 
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details. Recall that the bounce action 2ImS'[06] in Eq. (2.4) is the suppression exponent of 
false vacuum decay at zero energy, whereas the energy of the critical bubble gives the 
height of the potential barrier between the vacua. We also extracted the bounce residue Cb 
from the asymptotics of x) at —>■ —oo, see Eq. (2.1). We remind that the coupling 

constant (7 -C 1 scales out in the semiclassical calculations, cf. Eq. (1.3). 

We discretize equations (3.2) and introduce uniform Nt x 2Nx lattice with sites and xj 
covering the contour ABCD and space intervah^ {~L, L), L = 7; the spatial lattice spacing 
is Ax = L/Nx- We use the second-order hnite-difference approximation for the held equation 
and trade Fourier transform in Eq. (3.2c) for its discrete version. This turns the semiclassical 
boundary value problem into a set of Ntx2Nx nonlinear algebraic equations^^ for 0 ,^ = 0 s(fj, Xj) 
which are solved by the Newton-Raphson method [57]. Detailed description of our numerical 
technique will be presented elsewhere [58] , see Refs. [59] [36] for related works. In the subsequent 
Sections we will concentrate on numerical solutions at high energies and small multiplicities. 
We will need large W and because the typical frequencies Uk ^ E/N of these solutions are 
high. In particular, lattices Nt x = 3000 x 150 and 11000 x 4000 are required to reach 
acceptable numerical precision at E ~ Ecb and the highest energies, respectively. 

Our interest lies in the semiclassical results at high energies and relatively small N. The 
well-known saddle-point conhguration, however, is the bounce (j)b{t,x) which satishes the 
boundary value problem (3.2) at T = +cxd, 9 = 0 and E = N = 0. We continuously re¬ 
late it to the high-energy solutions. Namely, we hnd solntions at large T by iterative numerical 
method using the bounce as the zeroth-order approximation. After that we change T and 6 in 
small steps and obtain one slightly deformed solution (f)s(t,x) of Eqs. (3.2) at each step. We 
compute the energy E and initial particle number N for every 0s(f, x) by Eqs. (3.3). An exem¬ 
plary set of numerical solutions is marked by points (E, N) in the central plot of Fig. 5. It was 
obtained by decreasing T aX 9 = 0 (red/dark points), then increasing 9 to 9 = 0.4, and hnally 
decreasing T (green/light points) nntil the solution with T = 0 and 9 = 0.4 is reached (squared 
point). This procedure brings us to the high-energy region of interest. The other solutions in 
this region (not shown in Fig. 5) are obtained in similar manner, by changing (T, 9) in small 
steps and numerically solving the semiclassical equations (3.2). 

Solutions 0s = x) with T = 0 and arbitrary 9 are called “real-time instantons.” A 
representative solution of this type corresponds to the squared point in Fig. 5. In the next 
Section we will hnd that the real-time instantons are radically different from the solutions with 
T 7 ^ 0. We will develop an adequate high-energy description of collision-induced tunneling 
based on their properties. 

At low energies our semiclassical solutions reproduce the results of the previous Section. 
Indeed, Fig. 4c compares the semiclassical exponent (3.4) at 9 = 0 (points) with the exponent 
in Eq. (2.9) (line), where the numerical values in Fig. 4b are used. The two graphs coincide 
at E Ecb despite the fact that the perturbative results are obtained at N = 2 whereas 
the semiclassical solutions involve large multiplicities N ~ 1 / 77 ^, see Fig. 5. This is because 

Larger intervals are used at low energies due to larger sizes of the respective solutions. 

^'^The solutions are P-symmetric, 0s (h x) = 0s(i, —x), so we use only a half of the lattice with Xj > 0. 
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Figure 5. Central plot: (E, N) plane of initial data. Parameters of the numerical solutions are 
indicated by red (dark gray), purple and green (light gray) points. Circles with letters represent 
solutions from the insets (a)—(d). The squared point corresponds to the real-time instanton with 
T = 0 and 9 = 0.4. Insets (a) — (d): Three-dimensional plots of the numerical solutions Iie(ps{t,x), 
where the abscisses Ret — Imt parametrize the contours in the complex time plane. Color shows 
aig{4>s — 1). Euclidean parts of the solutions are marked by the square brackets with arrows. Above 
each 3D plot we draw the respective contour in the complex time plane, where the zoomed circular 
area with magnification factor displays the vicinity of the nearest singularity (crossed circle attached 
to the branch cut). The lattice size of all solutions is Nt x Nx = 3000 x 150. 
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Fn{E) is independent^® of at £' < Ecb [46], so we do not have to continue it to g'^N —)■ 0. 
Importantly, the perturbative graph in Fig. 4c becomes negative at ~ 0-8Ecb indicating 
apparent violation of unitarity. At these energies the perturbative expansion of Sec. 2 is not 
reliable unlike the semiclassical method of this Section. 

Our numerical results show a dramatic change in the form of the saddle-point solutions 
at ~ Ecb- The low-energy solutions in Figs. 5a,b resemble the bounce: they contain long 
Euclidean parts and describe creation of true vacuum bubbles (the latter are shown by red 
(dark) in the hgures). Waves in the left parts of the plots represent initial particles. In 
contrast, aX E > Ecb the initial waves are sharper and the bubbles are smaller, see Figs. 5c, d. 
Another property of high-energy solutions is small durations of their Euclidean evolutions and, 
nevertheless, nonzero values of the suppression exponents^® (3-4) thanks to complex-valued 
(j)s{t,x). In particular, the real-time instanton at T = 0 evolves entirely along t G M. We will 
argue that this feature guarantees stable perturbative expansion at high energies. 

One wonders how the durations of Euclidean evolutions can be of any physical meaning: 
they are not even explicitly specihed in Fig. 4a. Indeed, one expects that (j)s{t,x) are analytic 
functions of time and can be continued to any contours. But in fact, the solutions and, in 
particular, the bounce have branch-cut singularities starting at t which separate their 
time contours from the real time axis, see Fig. 4a. We compute positions of these singularities 
for the solutions (a)—(d) in Fig. 5. To this end we continue (j)s(t, x) to complex t and hnd 
the points t = where |E(0s)| = oo, see the zoomed areas in the respective t-planes. The 
singularities prevent us from continuing the low-energy solutions to the real time axis. We hnd, 
however, that Imt* decreases with energy and becomes negative for the solution (d) in Fig. 5. 
The latter solution can be considered in real time, as well as the real-time instanton at T = 0. 

We hnd the semiclassical solutions at g’^N = 2 (not shown in Fig. 5) and plot the respective 
suppression exponent E]\j{E) in Fig. Ic (thick line in the left part of the lower panel). The 
latter exponent monotonously decreases with energy until the point T = 0 and E = Ert{N) 
is reached^^. The region of few-particle initial states N cannot be directly addressed 

at the present-day computers. So, we extrapolate numerical results into this region. Since e“® 
analytically enters the semiclassical equations (3.2), the particle number N, action S'[0s], and 
the last term in Eq. (3.4) have regular Taylor expansions in this parameter. Moreover, e“® —)■ 0 
leads to Feynman boundary conditions at t —)■ —oo and therefore to g‘^N —)■ 0, see Eq. (3.2d). 
Substituting all Taylor expansions into Eq. (3.4) and using Eq. (3.5), one hnds that 

En{E) + g^9N + g^N = E{E) + 0{g^N^) . 

We ht numerical data for the left-hand side of this equality with function E{E)+d{E) ■g'^N'^ and 
obtain the suppression exponent E{E) (solid line in the left upper part of Fig. Ic). Numerical 

the thin-wall approximation which works well [45] at V{(j)+) = —0.4. 

^^Exponentially suppressed transitions at E > Ecb are called “dynamical tunneling” [60, 61] to distinguish 
from the potential tunneling at low energies. 

^^Data lines are not shown at high energies where the numerical errors are large. To decrease the com¬ 
putational cost, we find the solution at E = E^t at high precision (next Section) and fill the gap with thin 
interpolating line. 
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error of this procedure is expected to be smaller than 5%. Our results show that the two-particle 
exponent F{E) is also a decreasing function of energy. 

We hnish this Section by remarking that the above semiclassical solutions describe clas¬ 
sically forbidden transitions from the initial states with relatively small multiplicities N. In¬ 
creasing N, one reaches the states decaying classically (region “Fat = 0” in Fig. 5). Classical 
transitions in the model (3.6) were studied in Ref. [62], see Refs. [38-40] for the related work. 

4 Real-time instantons 

The real-time instantons, i.e. solutions of the saddle-point equations at T = 0 and arbitrary 6, 
are special in many respects (see one of them in Fig. 6a). At a given N they have the highest 
energies E = Ert{N) and smallest suppression exponents Fmin{N) = FN{Ert{N)) considered so 
far. Indeed, all solutions at smaller energies have positive parameter T: we obtained them by 
lowering T from T = -|-cx) to T = 0. Since OeFn = —2g‘^T, the exponent Fn{E) decreases with 
energy and reaches a local minimum at E = Ert{N). In Sec. 5 we will argue that this minimum 
is global by showing that the suppression exponent is energy-independent at F > Ert{N). 

The above property implies that each real-time instanton describes collision-induced tran¬ 
sitions from the states with hxed multiplicity N and arbitrary energies. Indeed, the probability 

dE , 

receives dominant contribution near the saddle point F = Ert{N) corresponding to the real¬ 
time instanton. One hnds ~ Q-Fmin(N)/g^^ where Fmin{N) is computed using x). 

In the previous Section we obtained the real-time instantons from the larger family of 
semiclassical solutions with arbitrary T and 9. One can compute (prt directly by solving the 
semiclassical equations (3.2a), (3.2b) with the initial condition 

Qk = e~%k , (4.1) 

cf. Eq. (3.2d). Then the parameter N, energy Ert{N) and minimal suppression exponent 
Fmin{N) are given by the standard expressions (3.3), (3.4). 

A convenient and important feature of the real-time instantons came as a surprising nu¬ 
merical fact in Fig. 5: they are dehned along the real time axis t G M. Technically, this prop¬ 
erty is related to the initial condition (3.2d) and its consequence (4.1). Indeed, the positive- 
and negative-frequency terms in the integrand of Eq. (3.2c) are of order and 

(ik = flfc ^ \vhere Tab = Iab is the height of the time contour. The integrals of these 

terms in Eq. (3.2c) converge faster and slower at higher Tabi respectively. At Tab = T the 
terms are of the same order because dk = bk due to the initial condition (3.2d). This “op¬ 
timal” contour lies right in the middle between the singularities of the solution: at somewhat 
higher or lower Tab one of the integrals in Eq. (3.2c) diverges signaling that the contour hits 
the singularity. We use Tab = T in numerical calculations and show the corresponding contour 
in Figs. 5a-d. At T —)■ 0 the optimal contour coincides with the real time axis. 
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Figure 6. (a) Real-time instanton (j)rt{t,x) and (b) its initial energy distributions at different Nx 
(numbers in boxes). In both figures 0 = 0.199, (F, N) ss (2.4Feb, 3.5/^^), Nt = 11000. Figure (a) 
uses colors of Fig. 5 and Nx = 2000. 


The above argument turns quantitative if we use the high-frequency asymptotics of the 
solution. Namely, consider the energy of modes with the wave number k in Eq. (3.3), 

> 9^^ = j dksk . (4.2) 

In Appendix B we demonstrate that any smooth solution has exponential asymptotics, 

^0 as k ^ +oo , (4.3) 


where T* is a parameter of the solution. Extracting ak and bl from Eq. (4.3) and initial 
condition (3.2d), one finds that the integral in Eq. (3.2c) converges at \Tab — T\ < T* and 
diverges otherwise. Thus, T* is the distance from the optimal contour to the nearest singularity 
of the solution, T* = |T — Imt*|. 

In Fig. 6b we plot energy distributions Sk for the real-time instanton with 6 = 0.199 at 
different Nx- At larger Nx and Uk m the graphs approach Eq. (4.3) with T* 0.013. This 
indicates that the continuum limit of our numerical solution is a smooth conhguration with 
singularities at hnite distances T* from the real time axis. 

Note that the latter fact is important because our conclusions about existence and proper¬ 
ties of the real-time instantons rely on numerical calculations. In Fig. 7 we study the continuum 
limit in more detail. If x) are smooth configurations, the numerical errors are expected to 
be polynomials in h = oc (Ax)^ and because we use the second-order hnite-difference 
methods. This is indeed the case: lattice values of E, A and Emin (points in Fig. 7) are not 
sensitive to A^ and well approximated by quadratic functions of 6 (lines). At large Nx all 
numerical errors are proportional to 6, see the inset in Fig. 7a. Once again we conclude that 
our lattice solutions with T = 0 have well-dehned continuum limits. 
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Figure 7. Continuum limits Nx —> +00 of the (a) real-time instanton energy Ert, (b) its initial 
particle number N and suppression exponent Emin] G = 0.199. We fit the data points with quadratic 
functions of (5 = N~‘^ (lines) in the range Nx = 800 4000. 




Figure 8. (a) Suppression exponent Emin{N) and (b) energy Ert{N) of the real-time instantons. 
Dashed lines in Fig. (b) are the lattice results at Nt = 11000 and different Nx (numbers near the 
lines). Thick solid lines are obtained by extrapolating results to Nx —)• + 00 . 


Numerical results for the minimal suppressions Fxnin{N) and respective energies Ert{N) are 
shown in Fig. 8. Dashed lines in Fig. 8b are the lattice results^®, while the solid lines represent 
continuum limits obtained by quadratic extrapolations in 5 = to 5 = 0, cf. Fig. 7. Points 
E = Ert in Fig. Ic are extracted from Fig. 8. In particular, we obtain results for the few-particle 
initial states g‘^N 1 by extrapolating Ert{N) to g'^N = 0 with linear function, and Ert{N) 
by the method described in the end of Sec. 3, see the dotted lines in Fig. 8. The accuracy 
of our result for F,.t(0) is better than 5%, while extrapolation of energy should be considered 

^®In Fig. 8a they are indistinguishable from the continuum limit. 
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Figure 9. (a) Numerical solution (j)s{t,x) at F ~ 2.QEcb > Ert{N). The other parameters are 
g'^N Ri 3.5, T Ri —2.8 • 10“"^, and = 1000. (b) Parameter T in the region E > Ert{N) for different 
Nx (numbers in boxes), (c) Energy distributions for the solution (a) and the real-time instanton with 
Nx = 1000 from Fig. 6b. In all figures 6 = 0.199 and Nt = 11000. 


as illustrative. In particular, we cannot completely exclude the possibility that E^t —t +oo at 
g'^N —)■ 0. Nevertheless, it is likely that the point E = Ert exists for the few-particle initial 
states, since it does for the multiparticle ones, cf. Refs. [25, 26]. 


5 Transitions at £” > Ert 

Numerical solutions sX E > Ert{N) look similar to the real-time instantons, cf. Figs. 9a and 6a. 
But in fact they are fundamentally different: we are going to demonstrate that they do not have 
continuum limits and lead to energy-independent suppression exponent E^^E) = Emin{,N). 

To begin with, we hnd that the lattice values of the param¬ 
eter T monotonously decrease with energy and become negative 
aX E > Ert{N), see Fig. 9b. Then Eq. (3.5) implies that the Eq 
exponent E]\[{E) reaches minimum sX E = Ert{N) and increases 
at higher energies. The last feature, however, is not expected 
in continuum models. Indeed, the sum in the dehnition (3.1) 
of Ei^{E) runs over all initial states with energy E and multi¬ 
plicity N. These include, in particular, the state where iV — 1 10 . Exclusive transi 

.. . tion at E > Ert(iV). 

particles perform transition at smaller energy Eq and one spec¬ 
tator particle carries the energy excess AE = E — Eq, see Fig. 10. The suppression exponent 
of the last process is En_i{Eq) ^ F;v(Eo), where correction of order g‘^9 ^ 1 is ignored, cf. 
Eq. (3.5). Since the inclusive transition is less suppressed than the exclusive one, we conclude 
that E]s[{E) < Eisf{Eo), i.e. E^^E) is a non-increasing function of energy. 

The above argument suggests that the negative values of T oc —OeE]^ in Fig. 9b are lattice 
artifacts and one should be careful with the continuum limit Ax —?• 0. Indeed, Fig. 9c shows that 
the high-frequency modes of solutions with E > Ert{N) are enhanced as compared to the case 
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of the real-time instanton and Eq. (4.3). To perform the quantitative comparison, we boldly 
assume that the semiclassical solutions at T < 0 have the form 0s(x) = (j)rt{x) + S(f){x), where 
6(j) consists of high-frequency modes evolving linearly in the real-time instanton background. 
Ignoring reaction of the latter on 6(j), we write, 

/ rik 

- - . (5.1) 

This solution is real at t —)• +oo because 5ck and 5c]. are the mutually conjugate constants. 
One expresses 6ck from the initial condition (3.2d): 5ck = — fli^*^)/(l ~ Ik)-, where 

Here and below we mark quantities related to the real-time instanton with 
“(rt)”, e.g. We obtain the energy distribution for the solution with T < 0, 




+ 5ck)ib 


(rt)* 

k 


+ 


(rt) sinh^(6'/2) 

^ sinh^(a;fcT + 6*/2) ’ 


(5.2) 


where /Att. Recall that Eq. (5.2) is based on a crude assumption that the 

solutions with E > Ert{N) differ from the real-time instantons only in linearly evolving high- 
frequency modes. Nevertheless, this scaling works well: in Fig. 9c the function (5.2) (points) 
coincides with the actual energy distribution (solid line). In Appendix B we derive Eq. (5.2) in 
somewhat different approach. 

Distributions (5.2) are qualitatively different at T > 0 and T < 0. In the former case 
the high-frequency asymptotics of Sk is consistent with Eq. (4.3) which implies smooth con¬ 
tinuum limit. At T < 0, however, the function (5.2) develops a nonintegrable singularity at 
cufc = cj* = —6/2T. Technically, this feature is related to the fact that the linear mode with 
frequency ca* satishes reality conditions at t —)■ ±cxd, see Eqs. (3.2b) and (3.2d). Solutions can 
accumulate macroscopic energy in this mode before its amplitude ~ joj^ becomes 

large and linear approximation (5.1) breaks down. Importantly, the latter energy tends to 
inhnity as T —)■ —0 or ca* —>■ -1 -cxd. 

Lattice solutions do not feel the singularity in Eq. (5.2) if the maximal lattice frequency 
^max = 2//S.X is below a;*. If we decrease /S.x with constant (T, 6*), the value of iXmax approaches 
cj* and the solution gains energy. Vice versa, at hxed energy one obtains smaller |T| at smaller 
Ax. If Eq. (5.2) is valid up to inhnitesimally small Ax^ parameter T approaches zero as Ax —)■ 0, 
so that ca* is kept outside of the lattice frequency range and the energy remains hnite. The 
respective semiclassical solutions arrive at the real-time instanton plus modes with inhnitely 
high frequency and vanishingly small amplitude carrying the energy excess E — Ert{N). These 
solutions do not have smooth continuum limits. 

The scaling property of our solutions is best summarized by dividing the total energy 
into “soft” and “hard” parts, i.e. E = Egoft + Ehard, involving modes with ojk < oja and 
uja < ojk < ojmax In Eq. (3.3). In Fig. 11 we plot EgoftiT) and EhardiT) for ua = 70m and 
several values of N^- Predictably, Egoft is not sensitive to while Ehard sharply depends 
on it according to Eq. (5.2) (lines in Fig. 11b). Since the data points in Fig. 11b are well 
approximated by Eq. (5.2), the value of T approaches zero as —)■ +C )0 at hxed Ehard, so that 
Egoft in Fig- Ha. 
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Figure 11. Energies of (a) soft and (b) hard modes of solutions at E > Ert{N) and different 
(numbers near the lines); 6 = 0.199. 

A remark is in order. Our numerical results show that at high enough E — Ert{N) the 
singularity ca* falls below Umax and interaction of high-frequency modes becomes important. 
We hnd, however, that the above qualitative picture holds even in this case. In particular, 
|T| decreases with at hxed total energy E, see Fig. 9b. We therefore expect that the 
respective continuum limit has the same properties: parameter T vanishes as ^ oo and 
lattice solutions approach the real-time instanton plus high-frequency modes. 

It is important to point out that the lattice solutions lacking smooth continuum limits at 
E > Ert{N), are still capable of describing tunneling transitions in that region. Namely, instead 
of discretizing the semiclassical equations (3.2) one can start from the lattice path integral for 
the cross section (3.1). The latter is an ordinary integral over W x 2Nx variables (pij. At ^ 0 
it can be evaluated in the saddle-point approximation, where the saddle-point values of (pij 
are the semiclassical lattice solutions. In this approach the existence of well-defined continuum 
limits of the latter solutions is irrelevant. The overall semiclassical results are reliable if the 
exponent En{E) has a limit at Ax —>■ 0. Then the lattice solutions point at the dominant 
mechanism of quantum transition. 

We have already argued that our lattice solutions with E > Ert{N) describe processes 
shown in Fig. 10: they tend to the real-time instanton with energy Ert{N) plus a few spectator 
particles in the form of high-frequency waves carrying the energy excess E — Ert{N). The 
respective suppression exponent is constant because T oc —OeE^ approaches zero as Ax —)■ 0. 
We finally conclude that En{E) is constant and equal to Emin{N) at energies above the threshold 
Ert{N), see Fig. Ic and cf. Refs. [24, 26]. 

6 Stability of perturbative expansion around the real-time instanton 

Since the inclusive cross section does not grow exponentially with energy at F' > E^t, one 
assumes that the terms of its perturbative expansion in around the real-time instanton are 
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also bounded as E ^ +oo. Then one can compute the collision-induced amplitudes and cross 
sections perturbatively at = 2 and arbitrary high energies. 

To test this property, we directly address the two-particle inclusive cross section of collision- 
induced false vacuum decay, 




( 6 . 1 ) 


where Uj —)■ Too, the initial state describes two particles with total momentum P = pi + P 2 = 
{E, 0) in the false vacuum Tq, each final state T f contains an expanding bubble of true vacuum, 
and we ignore prefactors. Let us check the possibility of extracting the cross section (6.1) from 
the perturbative Green’s functions in the backgrounds of the real-time instantons, so that the 
corrections to the perturbative calculations remain small as E ^ Too. 

In this Section we compute only the “factorized” contribution to (t{E) which was dominant 
and responsible for the exponential growth of the inclusive cross section in the Euclidean ap¬ 
proach of Sec. 2. We will see that this contribution is exponentially suppressed if the real-time 
instanton is used as a background. A consistent perturbative expansion around (j)rt{x) will be 
developed in the next Section. 

We start from the Green’s function (1.4) between the false vacuum Tq and the dominant 
final state of the real-time instanton (j)rt{x)- In the leading order of the perturbation theory 
we substitute 0 with ^^(a; — Xq) in the integrand of Eq. (1.4) and find. 


Grt = 




/ 


(fxo Mxi 


Xq) • • • ^o) •> 


( 6 . 2 ) 


where^® 

= . (6.3) 

Note that the position xq of Grt is not fixed by the semiclassical equations (3.2). Nevertheless, 
the final-state wave functional ^^rtlGrt] depends on xq via Grt{x — xq). 

We use the LSZ reduction formula and turn Eq. (6.2) into the 2 —)■ n T 'Lrt amplitude. 
Gonsidering the initial particles, we trade two of ^’s in the integrand for their positive-frequency 
residues see Eq. (3.2c). The case of the final-state particles is more involved because 

at t —)■ Tc )0 the configuration Grt contains, apart from the outgoing waves, an interacting 
bubble. To handle this difficulty, we assume that the false vacuum decay occurs in a finite 
volume \x\ < L with periodic boundary conditions. Then at t —?• Tec the bubble fills all space 
and Grt{x) reduces to waves in the true vacuum, 

/ dh 

- ^ T as t ^ Too , (6.4) 

where k) are the on-shell momenta in the vacuum G+- The LSZ formula substitutes 

Grt{x — Xo) with its residue for every final-state particle. For clarity below we consider 

finite large L and do not study the limit L —)■ Too. 

^®Recall that S and d'o are independent of xq in the limit tij —^ Too, as we drop the terms oscillating with 
ti and tf. 
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We obtain, 


'^2^n+'^rt 


(Pxq 


h* b. 


Pi P2 
2-\-n 


9 


Ck,...Ck^Kt[4>rt] e-0-(A:i+-+A:.-pi-p2) _ 


(6.5) 


Here pj and ki are the momenta of the initial and hnal particles, the prime of A! reminds that 
only the factorized contribntion is considered. Unlike in Sec. 2, we do not explicitly integrate 
over xq because 'hrt depends on it in a nontrivial way. Indeed, since (j)rt{x) is real at t —)■ +cxd, 
its dominant hnal state is a coherent one [63] [46]: in the interaction representation 


= exp i f dk 1 l^+) > (6.6) 

where we introduced creation operators c|, in the true vacuum 4/+. Parameters Ck of 4/^ 
are precisely the hnal-state residues Ck in Eq. (6.4). One can extract dependence of the wave 
functional '^rt[(prt] on xq from its transformation properties c|, —)■ and —)■ 

under spacetime shifts x ^ x + Xq, where E+ = 2LU(0+) < 0 is the energy of the true vacuum. 

The amplitude (6.5) has almost factorized form. In Appendix A we integrate over the 
phase-space volume and obtain the inclusive cross section, 

a'{E) = \Art 6p,6p,T+[0,t]|2 [ d^X exp 1 zP ■ A - zP+A° + [ I , (6.7) 

J [ J J 

where P = Pi + P 2 = {E, 0). Equation (6.7) looks similar to Eq. (2.7) of Sec. 2, and yet it 
is entirely different. First, the initial-state factor {bp-^bp^l'^ oc decays exponentially with 

energy because higher momentum transfer from the initial particles to the soft background is 
less probable, see Eq. (4.3). Second, the hnal-state contribution is nontrivial because is not 
an eigenstate of energy, it absorbs diherent energies in diherent cases. Thus, the simple picture 
of converting all energy into the multiparticle hnal states with huge phase volume is lost. 

Since the exponent in Eq. (6.7) is large, the integral over A is evaluated in the saddle-point 
approximation. One hnds the extremum of the exponent A^ = (—2zT', 0) satisfying 

E = E++ f dk . (6.8) 

J Ang^ 

At E = Ert the solution is T' = 0 because the right-hand side of Eq. (6.8) coincides with 
the hnal energy of the real-time instanton, cf. Eq. (6.4). At P > Ert one obtains T' < 0. 
Substituting A* into Eq. (6.8), we hnd, 

where the suppression exponent is introduced in the second equality. 

Equation (6.8) implies that the suppression exponent of ct'(P) grows at P > Ert, 

dF' 

— = 2 /(T,-r)> 0 . ( 6 . 10 ) 


a'(P) = \ArtA>+[4>rtf 


exp 


< 2P(T' 


P) - 2E+T' + 
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In the next Section we will demonstrate that aX E = Ert the exponent E'{E) coincides with 
the trne semiclassical exponent E{E). The latter, however, is constant at high energies. Thus, 
the factorized contribution (6.9) is exponentially subdominant, E' > E, ai E > E^t- 

7 Perturbative method at high energies 

Before considering the dominant contribution, we remark that the perturbative approach of 
the previous Section is incomplete. Indeed, there is a large family of real-time instantons 
parametrized^'^ with N = Nq and position i/q. Every solution from this family has its own 
dominant hnal state which also depends on Nq and yo, cf. Eq. (6.6). In this Section we 
change notation N ^ No for the real-time instanton parameter, to distinguish it from the 
number of colliding particles in the process; from now on, the latter equals two. We stress 
that the values of Nq and yo characterize the background for the perturbative expansion. Their 
values are selected to achieve better convergence. 

Our receipt for perturbative evaluation of the collision-induced cross section at high energies 
and two initial particles is summarized as follows [29]. One starts from the (n+2)-point Green’s 
function (1.4) between the false vacuum Tq and the dominant hnal state of the real-time 
instanton — yo), cf. Eq. (6.6). One evaluates the path integral for the Green’s function 

perturbatively in the background of 0^: substitutes (j){x) = (j)l^°\x — Xq) + g6(j){x) and expands 
the integrand in g6(j). Integral over the would-be hat direction xo ^ yo remains in the Green’s 
function, cf. Eq. (6.2). One hnally extracts the perturbative amplitudes A2^n+^rt from the 
LSZ formula, turns them into the inclusive cross section with the machinery of Appendix A, 
and takes the limit No —)■ 0. The result is a perturbative expansion for u{E) which, as we are 
going to argue, is applicable at arbitrary high energies. 

Let us explain the role of the auxiliary parameter No characterizing the background solu¬ 
tion. On the one hand, recall that the real-time instanton describes transition from the 
initial states with Nq particles: the limit N = Nq ^ ^ formally corresponds to the vacuum 
initial state. Indeed, at 0 —)■ -|-oo one simultaneously obtains Feynman initial condition —)■ 0 
for and Nq —)• 0, see Eqs. (3.2d) and (3.3). Thus, the real-time instanton with Aq = 0 
and Xo = yo is the formal saddle-point conhguration^^ for the path integral (1.4), and the per¬ 
turbative series around it constitute the ordinary saddle-point expansion. On the other hand, 
the energy Ert{No) of the real-time instanton represents the extremum, SeEn = 0, and stays 
nonzero in the limit No —)■ 0, cf. Fig. 8b. This means that is singular at Aq = 0 because 
its typical frequencies Uk ~ E/No are inhnite. We therefore develop perturbative expansion 
around the smooth conhgurations with Aq > 0 and send Aq —)■ 0 in the end of calculations. 

^^Recall that we work in the finite-L box which explicitly breaks Lorentz symmetry. Boosted real-time 
instantons should be taken into account in the infinite-volume limit. 

^^Recall that (j)rt extremizes the classical action and, by construction of 'krt, at xq = yo serves as the saddle- 
point configuration for the integral with the final state in Eq. (1.4). Note also that the two initial particles of 
the process are represented by the two Abactors in the integrand which do not change the initial saddle-point 
conditions Ofc = 0. The latter are satisfied by at Aq —)■ 0. 


21 



I 

a; 


1 


0.5 


0 




g^No g'^No 


Figure 12. Quantities (a) e“® and (b) T* characterizing the real-time instanton as functions of its 
parameter Nq. Numerical data at Nt x Nx = 11000 x 3000 (solid lines) are fitted with quadratic 
polynomials in Nq (dotted lines). 


In Fig. 12 we plot parameter e“® of the real-time instanton and distance T* to its closest 
singularity as functions of Nq. Numerical data (solid lines) are well fitted by quadratic polyno¬ 
mials with zeros at A^o = 0 (dashed lines). The graphs support our expectation that as Nq —)■ 0, 
the real-time instanton tends to a singular conhguration with vacuum initial conditions. 

Expanding the integrand of Eq. (1.4) in g6(j), we obtain Feynman rules involving points, 
propagators and vertices, 



which explicitly depend on Nq and xq. At nonzero No or xo ^ yo one also obtains the tadpoles, 
i.e. terms in action proportional to 6(j), which characterize deviation of the background solution 
from the true saddle-point conhguration. The tadpoles coming from the terms at t —>■ — cxd 
vanish as Nq —)■ 0, and we do not consider them in what follows. The hnal-state tadpoles 
are related to the fact that is the dominant hnal state for the conhguration (j)rt{x — yo) 
which is diherent from our background (j)rt{x — xq). We will discuss them in the end of this 
Section. Since the elements in Eq. (7.1) explicitly depend on x, energy and momentum are not 
conserved along the lines and in the vertices. Rather, the Feynman rules in momentum space 
involve structure functions depending on the momentum Q transferred to the background. We 
will see shortly that the real-time instanton consumes total energy ^ Ert from the initial 
particles. 

At zeroth order of the perturbative expansion one uses 0(x) = (l)lf°\x — xo) in Eq. (1.4) and 
obtains the diagram in Fig. 13a. Extracting the cross section and summing over the hnal states, 
one arrives at the contribution (6.9) which is exponentially subdominant. Indeed, we argued 
in Eq. (6.10) that the suppression exponent F'{E) of this contribution grows with energy at 
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Figure 13. Diagrams for perturbative expansion around the real-time instanton. 

E > Ert- Besides, at E = Ert and t 0 we have T' = T* = 0 and therefore 

/ dk \ri 

= E{Ert) + O(iVo) . (7.2) 

AtiuY’ 

Here we used Eqs. (6.3) and (3.4), vacuum wave functionals 

'ho[ 0 ] = exp ^;^ 0 (A;) 0 (-/c)| 

and 4/+[0], spatial Fourier transform (j){k) of conhguration 0, and representations (3.2c), (6.4) 
of the real-time instanton^^. Since the dominant exponent E{E) is constant at high energies, 
we repeat that E'{E) > E{E) at E > Ert, i-e. the factorized diagram in Fig. 13a is negligible. 

The dominant diagram at F > Ert is shown in Fig. 13b. It describes propagation of two 
initial particles which transfer momentum Q = qi + q 2 to the background. The respective 
transition amplitude is^^, 

= Artnin -1) I Kt[4>rt] ^ (7 3 ) 

cf. Fq. (6.5). Here qt = pt — ki and Dp, is the double residue of the propagator {6(j) S(f))rt- 

In Appendix A we convert the amplitude into the inclusive cross section, 

r m 2 f dk\dk 2 ,2 

O'(D) = \Art / (ip) ppiy \Dpi,ki Dp2^k2 + -Opi,fc 2 -Dp 2 ,fci I 

J SZTI 

X [ d^X exp < iQ ■ X — iE+X^ + ( ^ A) 

J { J J 

where X = yo — xq and we omit trivial prefactors. The hrst line in this expression is the naive 
Feynman diagram in the background of the real-time instanton. The factor in the second line 

^^Recall also that Fi^{E) = F{E) — g'^BN + 0[N) as we argued in the end of Sec. 3. 

^^One finds that as L — 00 , the energy of the true vacuum becomes infinite and therefore T' —0. This 
means that in the true infinite-volume limit (not considered here) the contribution (6.9) becomes comparable 
to the dominant one. 

^'^Symmetrization with respect to permutations of fci, ..., A:„ is assumed. 
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is related to the on-shell hnal state of the background process. It is the same as in Eq. (6.7). In 
Sec. 6 we demonstrated that this factor is sharply peaked around = {Ert, 0) with fluctuations 
of order AQ gErf 

Evaluating the integrals over Q and A in the saddle-point approximation, we get with 
exponential precision, 

a{E) = at E> Ert , 

see Eqs. (6.9) and (7.2). We thus obtained constant suppression exponent E = E{Ert) at 
E > Ert which was deduced in Sec. 5 on the basis of sophisticated numerical analysis. Note 
that our perturbative method can be applied for calculating the prefactor: one just has to 
estimate the energy-independent saddle-point determinant in the Green’s function (1.4) and 
collect few simple prefactors in the above calculations. 

The hrst perturbative correction to the dominant contribution is shown in Fig. 13c. It 
involves the same relative factor g'^n^ as in the low-energy calculation of Sec. 2. This time, 
however, the number n of the point-like factors (f)rt representing the particles added to the 
background, is relatively small and does not grow with the collision energy E. Indeed, any 
energy transfer to the background above Q = Ert is cut off by the hnal-state factor. This 
means that n ~ |Q — Ertl/m ~ g~^, and the diagram in Fig. 13c is of the same order as the 
dominant one. Importantly, it does not grow with energy. Resummation of these disconnected 
diagrams [51, 52] leaves us with corrections involving vertices which are suppressed by the true 
expansion parameter g. 

Let us hnally discuss the tadpoles. We saw that the parameter uq of ^rt does not coincide 
with the position xq of the background solution: the integral over the difference \ = uq — xq 
enters Eq. (7.4). Thus, (j)rt{x—xo) is not the true saddle-point conhguration of the integral (1.4), 
perturbative expansion around it starts from the linear term in 6(j). We obtain the tadpole 

(f 

which should be integrated over q with the hnal-state residue of the propagator. Apart from 
the additional integration, this tadpole is similar to the hnal on-shell particle in the amplitude. 
One may think that the dominant contribution includes diagrams with the tadpoles attached 
to the hard propagators, e.g. Fig. 13d. However, all momentum pi = gi + g of the propagator 
with the tadpole is transferred to the real-time instanton. Indeed, Eq. (7.5) is proportional to 
the exponent which accumulates transferred momentum in the amplitude, cf. Eq. (7.3). 

Then the total transferred energy in the process in Fig. 13d is higher than Ert and the respective 
contribution to the cross section is exponentially small. 

We see that the tadpoles infest the hnal-state propagators in Fig. 13c and subdominant 
diagrams^^. Note, however, that the element (7.5) is 0{g^) at best. Indeed, at small q the 

^®One can check that the tedpoles do not change our conclusion about exponential suppression of the factorized 
contribution. 


,-iq-X 


1 ) 


(7.5) 
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bracket in Eq. (7.5) is proportional to g because X g ai E ^ Ert- At g ~ m/g the tadpole is 
suppressed by the factor which vanishes at high q because the energy of the real-time 

instanton is hnite. Thus, the tadpoles also do not break the perturbative expansion. 

We conclude that the perturbative expansion around the real-time instanton is reliable at 
arbitrarily high energies. 


8 Summary and discussion 

In this paper we studied collision-induced tunneling in held theory. We paid special attention 
to the case of two colliding particles with high total energy E. As a playground we considered 
induced false vacuum decay in (1 + 1) dimensions. We demonstrated that the suppression 
exponent En{E) of this process decreases with energy, reaches minimum E = Emin{N) at 
E = Ert{N) and remains constant at higher energies. 

Our methods rely on existence of the real-time instantons — a special class of semiclassical 
solutions describing inclusive collision-induced transitions from the initial states with N par¬ 
ticles and arbitrary energies. The minimal suppression Emin{N) and threshold energy Ert{N) 
are computed as functionals on these solutions. Real-time instantons were hrst observed in the 
toy model of Ref. [26]. Here we numerically obtained them in the case of (1 -|- l)-dimensional 
false vacuum decay. We expect that these solutions exist for other collision-induced tunneling 
processes. One can verify this expectation on a case-to-case basis by solving the respective 
semiclassical boundary value problem in a given model. 

Importantly, we argue on general grounds that the real-time instantons are complex solu¬ 
tions evolving in real time (hence the title). This property is very unusual for the semiclassical 
solutions related to exponentially suppressed transitions, it leads to far-reaching consequences. 
We hnd that scattering of high-energy quantum particles in the backgrounds of the real-time 
instantons resembles scattering in vacuum because energy exchange between the particles and 
the soft background occurs with exponentially small probability. This situation is radically 
different from that at low energies where the Euclidean semiclassical solutions recycle any 
additional energy into exponentially large probability factors. Starting from the real-time in¬ 
stantons, we develop a perturbative description of the two-particle collision-induced processes 
at high energies. We demonstrate that this description remains valid at arbitrary high energies. 
Our method shows that the suppression exponent E{E) = E 2 {E) is constant at E > Ert{‘2). 
The collision-induced transitions in this regime involve transfer hxed energy Ert from the two 
colliding particles to the soft background; the energy excess E — Ert remains in the initial par¬ 
ticles till the end of the process. Note that our perturbative methods can be easily generalized 
for calculating prefactors or exclusive cross sections. 

We conclude that the real-time instantons, if exist for a given collision-induced process, 
provide powerful perturbative framework and guarantee constant suppression exponent E{E) = 
Emin{‘^) of this process at energies above a certain threshold Ert{2). 
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A Multiparticle cross sections 


In this Appendix we evaluate inclusive cross sections for the amplitudes considered in the main 
body of the paper. We start from the factorized 2 —)■ n amplitude 


A 


_ Ao 

2—>-n ~ Cfcj . . . 

9 


(A.l) 


Here Aq and g are constants, fc* are momenta of particles in the hnal state. The cross section 
of inclusive transition to the n-particle hnal states is obtained by integrating over the phase 
space volume n„, 


an{P) = J \A2^n\^ dUn{P) 


|AoP 

f dki\ck^\ 

n\ j 

Ang'^Uk, 


dkn\ckS 


H- \- kn 


P), 


where we ignored the initial-state factor, introduced the total initial momentum and on-shell 
frequencies u‘l = k"^ + w?. We use Fourier representation of the ^-function and hnd. 


a„(P) = \^ [£\ [/(A)]"P^-" , where /(A) = [ ^^.2) 

n\ J J Airg^Uk 

If Cfc = Cfe does not depend on fc. 


/(A) = Po(W-A2 + ^eAO) . 

27rg^ 

Equation (A.2) shows that A is a typical Compton wavelength of the hnal particles, A ~ k~^. 
Summing up the n-particle contributions (A.2), we obtain Eq. (2.7) of Sec. 2. 

In Sec. 6 we consider the factorized amplitude (6.5) of the process 2 —)■ n + '^rt- To simplify 
summation over the hnal states in the inclusive cross section, we relate Eq. (6.5) to the ordinary 
2 —)■ n -|- m amplitudes. To this end we expand the exponent in the hnal state (6.6) of our 
amplitude and obtain. 


A' 


2—¥n-\-'^rt 


oo 

E 

m=0 


1 

ml 


dkn+i 


dk 


'n+m Ck 


Anguil^] 


A! 


■ 2 ^ 77+777 


(27r)V^)(P-P/) , (A.3) 


where Pj is the total momentum of the (n -|- m)-particle hnal state. Next, we compare the 
right-hand side of Eq. (6.5) with Eq. (A.3). Substituting 0 = 0^ into the wave functional Trf[0] 
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of a coherent state [28] with parameters h^e we hnd, 

{ r rj, p* p-ik-XQ 'I 

/ dk - \ . > e"*®+*° ^+[(j)rt] , (A.4) 

where we extracted dependence on xq as explained in Sec. 6, used the form (6.4) of took 
the limit t/ —)■ +C )0 -e"*^ and introduced wave functional 4/+ of the true vacuum. The hnal state 
'^rt in Sec. 6 has dk = Ck, cf. Eq. (6.6). Expanding the exponent in Eq. (A.4) and substituting 
the series into Eq. (6.5), we hnd, 


dkji+m ^ 

^n+m 

X "^Ck,... Ck„ {2Tif 5(2) (P - Pf) . (A.5) 

Here we evaluated the integral over xq and introduced = -A-rt ^,^^, 2 /9^- Comparing 

Eqs. (A.3) and (A.5), one hnds that the amplitudes ^ 2 ^.n+m have the form (A.l). We hnally 
obtain the inclusive cross section (6.7) by summing up the n-particle ones in Eq. (A.2). 

To process the dominant amplitude (7.3) of Sec. 7, we do the opposite to the above, i.e. 
combine the hnal states into the coherent states 4/^ with parameters dk = Cke^^'^° + dck, where 
cj^Qik-yo parameters of 4/^* and dck are arbitrary. We hnd. 


4' 


2—¥n-\-'^rt 


00 

E 

m=0 


m\ 


dkn+1 






E 

n=0 


nl 


dki 5c\^ 


= A 


rt 


d^xo T* 




dkn del 

_ A 

■^2— 


dki dk2 dc\^5c\^ 


F) 1 V) 1 p-^(^o-(Qi+92) 

■^pi,ki ■^P2,k2 ^ 


(A.6) 


In the hrst line of this expression we expanded in 6c using Eq. (6.6), in the second substituted 

Eq. (7.3) and introduced wave functional Eq. (A.4). 

Now, we evaluate the inclusive cross section 


17 ( f dk d' d* ^ 

a{E) = ^JVd' Vd* I “ y J ’ 

where E^^) is the spacetime volume. First, we note that the amplitude (A.7) depends on 
the arbitrary parameter yo. This is the freedom of choosing the background for expansion, it 
disappears after resummation of perturbative series. We hx the freedom requiring 


Vo = A , y'o = XQ , 


(A.8) 


where y'^ and x'q come from the complex conjugate amplitude (A.6) in the integral for the 
inclusive cross section. Second, the parameters 5c^ in the prefactor can be obtained by varying 
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the exponent exp{J dk akScl/{A7rg'^uj^^)} with respect to a at a = 0. The integral in Eq. (A.7) 
therefore can be evaluated using the functional 


n[Q;,a'*] = J V6c'V6c* '^d’[(j)rt] exp 


dk 


,(+) 




+ [(prt\\ exp 


dk 


A'ng'^oj^ 


(+) 


(afcCTfc + Icfcpe"*'^'-^) 


{akdcl + a'^dd^- d'^dl) 


(A.9) 


where we denoted d'^ = Cfce*^'^o + 5c'^, substituted Eq. (A.4), evaluated the integrals over 5c*, 
5c' and denoted \ = Hq — Xq = x'q — Xq. We see that variations over and a'^ give contraction 
rule for the hnal-state variables 5ck^ 5c\_^ = 5{ki — k 2 )- 

Substituting Eqs. (A.6) and (A.9) into Eq. (A.7), one obtains the inclusive cross section 
(7.4) containing the integral over \ = x'q — xq; integration with respect to (xq + x'q)/2 gives the 
two-volume which is canceled in Eq. (A.7). 


B High-frequency tail of the semiclassical solution 

Let us evaluate high-frequency asymptotics of the saddle-point solution (j)s{x). To this end we 
represent (j)s{x) = (t)o{x) + 5(j){x) as a sum of soft nonlinear background 0o and high-frequency 
part dcf) 00 evolving on top of it. Functions 0o and 50 contain modes with k < A and k > A, 
respectively. One rewrites the held equation as 


□50(a;) = J{x) , 


(B.l) 


where J is a contribution of 0o and 50 is ignored in potential terms. We solve Eq. (B.l), 


50(x) 


d^k J{k)e 

(27r)2 


k>A 


(B.2) 


using the two-dimensional Fourier image of the source J{k) and arbitrary on-shell waves in 
the second term with /c^ = {uk, k) and Uk = \k\. The solution (B.2) is real as t ^ +cxd in 
accordance with Eq. (3.2b). In the inhnite past, i.e. t ^ iT — oo, Eq. (B.2) takes the form 
(3.2c) with 


ak = Ck + Jk , K = <^k + Jk^ and J^ = ±iJ{Tk) 
We hnally solve the initial condition (3.2d) and obtain, 

^k{j^y - jy 


fc0=l 


(B.3) 




— 


1 - 7fc 


^ — „-2a;fc'r-e 

7fc = e 


(B.4) 
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Note that Eq. (B.l) and its solution (B.2), (B.4) are valid only at 
k > A ^ m. 


/\ 


© 


One notices that the sources are exponentially sensitive to J 
k. Indeed, the t-contours in the Fourier transforms 





k 



(B.5) Jfc 



can be deformed into the upper and lower parts of the complex 



time plane until they hit the singularities and t' of the solu- ^ ^ ^ ; r t+ 

Figure 14. Contours for 

tion"^°, see Fig. 14. At large k the integrals receive the dominant 
contribution near the singularities, and 



(B.6) 


where Jq are constants. 

We hnally compute the energy of the initial particles with momentum k in Eq. (3.3) using 


Eqs. (B.3), (B.4), 



(B.7) 


47r lQnsmb?{ujkT + 6/2) ’ 


Estimating the largest term in the nominator by Eq. (B.6), one obtains the high-frequency 
asymptotic (4.3) with T* = T + min(—Imf*, Imf(). Equation (B.7) obeys^^ the rescaling 
property (5.2) of Sec. 5. 
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